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We present a multiscale quantum-defect theory (QDT) for two identical atoms in a symmetric 
harmonic trap that combines the quantum-defect theory for the van der Waals interaction [B. Gao, 
Phys. Rev. A 64, 010701(R) (2001)] at short distances with a quantum-defect theory for the 
harmonic trapping potential at large distances. The theory provides a systematic understanding of 
two atoms in a trap, from deeply bound molecular states and states of different partial waves, to 
highly excited trap states. It shows, e.g., that a strong p wave pairing can lead to a lower energy 
state around the threshold than a s wave pairing. 
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I. INTRODUCTION 

Two interacting particles under confinement, described 
generally by a Hamiltonian 



H = -- 
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where V± and V 2 are the confining potentials and v(r) 
is the interaction between particles, represents a funda- 
mental class of problems in physics. One famous ex- 
ample is the helium-atom problem that has played an 
important role in our understanding of electron corre- 
lation in atomic physics (see, e.g., Ref. Similarly, 
the problem of two atoms in a harmonic trap, which 
has attracted considerable recent attention (see, e.g., 
Refs. @, i 1, i, 0, 0, 1, & M HO, M El ) , is the key to our 
understanding of atomic correlation in a trapped many- 
atom quantum system. Such correlation differs qualita- 
tively from the electron correlation because atoms attract 
each other at large distances and can form bound states. 

Existing theories of two atoms in a trap have relied 
mostly upon the pseudopotential model of atomic inter- 
action [14[ , and its generalizations [7| . While such models 
can work well in describing how the trap states, espe- 
cially the lowest few, are affected by atomic interaction, 
they generally fail in describing how a molecular state is 
affected by trapping, with the only exception being the 
least bound molecular state with a very large scattering 
length. Furthermore, such theories do not adequately 
address nonzero partial waves, for which naive general- 
izations of the "shape-independent" approximation [2J, 
using, e.g., the effect range theory (ERT) (l5| . would 
generally lead to incorrect results. 

We present here a multiscale (2 length scales, to be 
exact) quantum-defect theory (QDT) for two identical 
atoms in a symmetric harmonic trap. R is a completely 
general theory that works for different partial waves, and 



from deeply bound molecular states to highly excited 
trap states. In Sec. UH we expand our tool box of QDT 
for different long range potentials [3, [M E3, H El by 
presenting a QDT for a symmetric harmonic potential. It 
is independently useful beyond the scope of two atoms in 
a trap. For example, it may be used to treat two nucleons 
outside of a closed shell [2fj| . 

In Sec. IIII1 this theory is combined with the angular- 
momentum-insensitive quantum-defect theory (AQDT) 
for the van der Waals interaction [2l| to formulate a two- 
scale QDT that provides a systematic understanding of 
two identical atoms in a symmetric harmonic trap. Re- 
sults and discussions are presented in Sec. lIVi including a 
discussion of the limitations of "shape-independent" ap- 
proximations, and a universal spectrum for two atoms in 
a trap at the van der Waals length scale [2l|, [H, [H, [24j 
that shows, e.g., that a strong p wave pairing can lead 
to a lower energy state around the threshold than a s 
wave pairing. We will also show that two atoms in a 
trap has long-range correlation that becomes important 
for large scattering lengths, a result that has proven to 
be the key for generalizing the variational Monte Carlo 
(VMC) studies of few atoms in a trap to the regime of 
strong coupling [24|. Conclusions are given in Sec. [V] 



II. QUANTUM-DEFECT THEORY FOR A 
SYMMETRIC HARMONIC POTENTIAL 

The goal of a QDT for a symmetric harmonic potential 
is to provide a systematic understanding to a class of 
problems described by the radial Schrodinger equation 
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Unlike the standard textbook solution which requires 
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to as the "pure" harmonic oscillator, the QDT formula- 
tion is applicable to any V{r) that is asymptotically a 
harmonic oscillator, but may differ from it at short dis- 
tances in an arbitrary fashion. 

As in any QDT formulation [j], [l(| [l7j], we start by 
defining a pair of reference functions that are two lin- 
early independent solutions for a symmetric harmonic 
potential 
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The solutions can be easily found [I 



v(r) = . (4) 
and we will take 
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Here e = e/fko is a scaled energy, x — r/(3ho is a ra- 
dius scaled by (3ho = (^/a*^) 1 ^ 2 , which is a length scale 
associated with the harmonic potential. M is the con- 
fluent hypergeometric function [251, b = (I + 3/2 — e)/2, 
and c = I + 3/2. In this definition, fj£°* is regular at 

the origin, and g^°^ is irregular. They will be called the 
regular solution and the irregular solution, respectively. 
They are also chosen such that their Wronskian is given 
by 
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With this definition of reference functions, the wave 
function u e i for any potential that is asymptotically a 
harmonic oscillator can be written, at sufficiently large 
distances, as 

u d (r) = A4f^°\x) - KM(e,l)g$°\x)] . (8) 

This defines the K matrix, K^ ho \ for a symmetric har- 
monic potential, with its value being generally deter- 
mined by matching Eq. ([8]) to the short-range solution. 

Making use of the large r asymptotic behaviors of /C 10 ) 
and g( ho \ as given by [25| 
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and enforcing the boundary condition, u e i(r) — ► 0, at 
infinity, we obtain the following equation that gives the 
the energy spectrum of Eq. @ as the crossings points of 
two functions: 



xl M (e)=^°)(e,0 



(11) 




FIG. 1: The x\ H °\e) function for the s wave. The energy 
spectrum for any potential that is asymptotically a symmet- 
ric harmonic oscillator is given by the crossing points of this 
function with a short-range K matrix K {ho) (e,l) defined by 
Eq. ©. The special case of K (ho) (e,l) = for all energies 
corresponds to a pure harmonic oscillator. 




FIG. 2: The same as Figure 1 except that it is for the p wave. 



Here 



X f OJ (e) = (-l) i+1 [r(/ + 3/2)] 
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(12) 

is a universal function of the scaled energy e that is deter- 
mined solely by the harmonic potential, and K^ ho ^ is the 
K matrix that encapsulates all the short-range physics. 

Plots of xt°^ ( e ) l° r t ne s an< l P waves are shown in Figs. [I] 
and [21 respectively. 

Equations (jTTJ) and (fT2|) give a rigorous formulation 
of energy spectrum for any potential that is asymptoti- 
cally a harmonic oscillator. It can be implemented nu- 
merically, or used as a basis for approximate analytic 
solutions. The case of a pure symmetric harmonic os- 
cillator is included as a special case, corresponding to 
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X\ h °)[e,l) = for all e, with a well-known energy spec- 
trum of e = hu(2j + I + 3/2), where j = 0, 1, 2, . . .. 
From Eq. (jTTJ) , it is clear that the key to finding the 

energy spectrum is to find the K matrix K^ ho \ as xt°^ 
is already known analytically. For this purpose, we note 
that the reference functions /C 10 ) and g( ho ^ have the fol- 
lowing small r asymptotic behaviors that can be derived 
from a proper expansion of the confluent hypergeometric 
function 12a. For e > 0, we have 



trap reduces to the solution of Eq. ((2|) with 
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where k — (— 2(ie/h 2 ) 1 / 2 . For interactions V(r) that de- 
viate from the harmonic potential only in the region of 
t <C Pho, these behaviors greatly facilitate the match- 
ing to the short-range solution, from which K^ ho ^ can 
be determined. This is illustrated in our treatment of 
two identical atoms in a symmetric harmonic trap, to be 
presented in the next section. 



III. TWO IDENTICAL ATOMS IN A 
SYMMETRIC HARMONIC TRAP 

Two interacting atoms in symmetric harmonic trap are 
described by a Hamiltonian 



V(r) 
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where represents the radius inside which the inter- 
actions of shorter range than fie, such as the — Cs/r 8 
correction, would come into play. 

Since the V(r) characterized by Eq. (fi"9|) is asymptot- 
ically a harmonic oscillator, it is amenable for the QDT 
treatment of Sec. HU In particular, the solution for the 
energy spectrum reduces to finding the K matrix K^ ho ^ 
for the class of problems defined by Eqs. (fl~8|) and (fT9j) . 
This will be accomplished here by taking advantage of 
the disparate length scales in the system. 

For r > r , the potential V(r) as given by Eq. (flU)) 
has two length scales. In addition to (3q that is associ- 
ated with the van der Waals interaction, the harmonic 
trapping potential has a length scale that can be taken 
either as the (3ho defined earlier, or as a/ lo = (Ti/muj) 1 / 2 = 
fiho/V?*- We will use both interchangeably, but will em- 
phasize aho for the sake of easier comparison with other 
results. 

For atoms in a typical magnetic or optical trap, a^o is 
of the order of a micron, which is much greater than (3q 
that is of the order of 100 au or about 50 nm. Under 
this condition of (3q <C a/ lo , which we will call the limit 
of weak confinement, the van der Waals and the trap- 
ping potentials operate on distinctive lengths scales. In 
the region of r ^ fie, the van der Waals interaction is 
negligible, and we have 



u u {r)=A[f^-K^{e,l)g^] 
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In the region of ro < r <C a/ lo , the harmonic potential is 
negligible, and the interaction is dominated by the van 
der Waals interaction. Here the wave function can be 
written as I2l| 
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where v(r) represents the interaction between them. To 
be specific, we restrict ourselves here to a class of prob- 
lems for which the atomic interaction is characterized, 
at large distances, by an attractive 1/r 6 van der Waals 
potential 



v(r) 



-C 6 /r 



(18) 



which has an associated length scale of f3g — 

(2^/n 2 ) 1 ^ m. 

For two atoms having the same trapping frequencies, 
namely io\ = u)2 = which include of course the case of 
two identical atoms of interest here, the center-of-mass 
motion and the relative motion are separable, and the 
solution of two identical atoms in a symmetric harmonic 
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where /° j and g° e j are the reference functions for the 
van der Waals potential, and K c is the corresponding 
short-range K matrix [lj, H2, [H| . 

For weak confinement defined by fi@ -C aho, there ex- 
ists an intermediate region /3g <C r <C a/ lo in which either, 
or both, the van der Waals potential and the trapping po- 
tential can be ignored. In this region, the wave function 
can be written either in the form of the inner solution as 
given by Eq. pTjl , or in the form of the outer solution as 
given by Eq. (|2TJ|) , and they must agree with each other. 

Since r/ai w <C 1 in the intermediate region, f^ 1 "^ and 

<7g2 are given by Eqs. (fll?|) - (fll)|) . In the same region, 
r/Pe ^> 1 and the reference functions for the van der 
Waals potential are given for e > by [H, [26| 
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and for e < by pj, 

flf = i[M/f ) +2(-l)^;f](r//3 6 ) 1 / 2/;+1/2 ( Kr) 

+ l\wf> 2(-iywff](r/(3 6 )^I_ l _ 1/2 ( K m) 
9lf = i[^ 6) +2(-l)'<f](r//3 6 ) 1 /2 /;+1/2 ( Kr) 

+ \[wf ] 2(-l) , W fl fJ](r/Ai) 1/2 J- I -i/ a (/c«5) 
I 

< 6) ( es ) 
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Here M e ; = G d {-v)/G d {v), with z/, X ci , F e ;, and G dl 
all of which are functions of the scaled energy e s , being 
defined in Ref. (l8| . 

Comparing, in the intermediate region, the inner so- 
lution given by Eqs. (EH) and (|2"2 ]) -([2"5 ]) with the outer- 
solution given by Eqs. (f20f and (fT3|) -(|16 | . we obtain, for 
e > 0, 

i 



Here the Z c ^> and W c ^ matrices describe the propa- 
gation of a wave function in a —Cg/r 6 type of potential 
from small to large distances, and vice versa Hi], [2tJ ■ 
Their elements of are all universal functions of a scaled 
energy e s = e/s_e, where se = (h 2 /2fj,)(l/ (3q) 2 is the en- 
ergy scale associated with the van der Waals interaction. 
Explicit expressions for the elements of the Z c ^ can be 
found in Ref. [13] ■ The W c ^ matrix, which is related to 
the W matrix defined in Ref. [18[ by a linear transforma- 
tion [HI], is given by 



(26) 
(27) 
(28) 
(29) 



where 



tan^e) = -{ Z f? - Z «*> K<){Z c f f - Z c g f , 

is the physical K matrix for atomic scattering in free 
space as given in AQDT [2l|. For e < 0, we obtain 
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where X; (es) = Wf_ /Wg_ is the %-function that de- 
termines the molecular spectrum in the absence of trap- 
ping [m. 

Equations [[IT]) - ([12]) and ([3"0 |) -([3"2 ]) give a complete de- 
scription of the energy spectrum for two idential atoms in 
a symmetric harmonic trap, from deeply-bound molecu- 
lar states to highly excited trap states. The only assump- 
tion in the theory is the assumption of weak confinement 
as specified by fte/aho 1, which is well satisfied under 
all existing experimental conditions. 



Other than the two energy scaling parameters fuv and 
Se, the only parameter in the theory is K c (e,l), which 
characterizes the interactions of shorter range than (3$. 
It can be replaced by other equivalent short-range pa- 
rameters such as t he q uantum-defect /i c (e, I) [29( or the 
i^;°(e,0 parameter [la [n|. 



5 



IV. RESULTS AND DISCUSSIONS 



A. Universal spectrum at the length scale of /3g 



The results given by Eqs. (HU, (JT2J, ([50 ]) -([52 ]) can be 
written in different forms that are convenient for differ- 
ent purposes. For conceptual understanding, it is best 
to rearrange them so that the entire energy spectrum is 
given by the solutions of a single equation 



xl h °' 6) (e,p 6 /a ho )=K c (e,l), 



(33) 



where 



(ho,6) = Zf:\e s )-Z^{e s )^\e) 



for e > 0, 
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for e < 0, and we have defined 
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The xi h0 '^ function is a universal function that is ap- 
plicable to any two identical atoms in a symmetric har- 
monic trap, provided they interact via the —Cq/t 5 type 
of van der Waals potential at large interatomic separa- 
tions. The strengths of interactions, as characterized by 
Cq and uj, play a role only through energy scaling pa- 
rameters se and %uj. Specifically, the Xt function is 
made up terms that depend on energy through two dif- 
ferent scaled energies, e s = e/sE and e = e/hco that are 
related by e s = (/3e/aho) 2 ^- In other words, it is made up 
functions that vary on two distinctive energy scales: the 
anc | the a^ ho ' functions that varies on the scale of 
Tioj, and the Z c ^ and W c ^ matrix elements that vary 
on a scale of se ■ 

The solutions, more precisely the inverse, of Eq. (|33|) 
can be written as 



hf\K c ,j3 6 /a ho ) 



(38) 



where, similar to X; i the flf are universal functions 
that are uniquely determined by the exponent of the van 
der Waals interaction (n = 6), and the exponent of the 
trapping potential (2 for the harmonic trap). Instead of 
the parameter K c , the same universal functions can also 
be expressed in terms of equivalent parameters such as 



the quantum-defect \i c [29| , or the Kf parameter [3(| , all 
of which are well defined for all /. 

When the energy and angular momentum dependence 
of the short-range parameter K c is included, Eqs. (|3"3"]) 
and (|38|) are exact, and applicable to arbitrary energy 
and partial waves. Ignoring the energy and the I de- 
pendence of K c , which in the case of a single channel 
is due entirely to interactions of shorter range than /3g 
[H|, the solutions of Eq. ([33]), namely Eq. (38]) with 
K c = K c (e — , / = 0), or its variations, give what we call 
the universal spectrum at length scale (3§ [2l|, ]22j, ]24j . It is 
followed by all two-atom systems in a trap with —C'e/r 6 
type of long-range interaction, over a range of energies 
that is hundreds of se around the threshold [H], EH, 
which far exceeds all energies of interest in cold-atom 
physics. (As an example, se = 9.331 x 10 -4 K for 23 Na.) 
Other than the energy scaling parameter se that is de- 
termined by the Cq coefficient and the atomic mass, all 
energy levels in this energy range, including states of dif- 
ferent I, are determined, in the case of single channel, 
by two parameters, K c = K c (e = 0, 1 = 0) and (3e/aho- 
This is an example of universal spectrum at the second 
longest length scale in the system (24|, as opposed to 
the universal spectrum at the longest length scale, which 
would have required only a single parameter (2lj |. 

Depending on the physics of interest, the universal 
spectrum at length scale /3@ can also be expressed in terms 
of other parameters. In particular, it can be expressed in 
terms of the s wave scattering length, ao, as 

e = £li(a /a ho ,P 6 /aho) , (39) 
since ao can be related to K c = K c (e = .1 = 0) by 



■ ■2b 



T(l-b) 



T(l + b) 



K c (0,0) + tan(7r&/2) 
^ c (0,0)-tan(7r6/2) 
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where b = l/(n — 2) with n = 6 [29L 1321 ] . Similar repre- 
sentations of universal spectrum can also be defined more 
generally for N atoms (N > 2) in a trap [22l. [23l [24j . 

Figure [3] illustrates the universal s wave spectrum at 
length scale (3q for two identical atoms in a symmetric 
harmonic trap. It uses the representation of Eq. (|39|) 
to facilitate comparison with the shape-independent ap- 
proximation (see Ref. 0] and Sec. HVB]) . The f^ z=0 ) func- 
tions of two variables a^/aho and Pe/ a ho are plotted here 
as functions of a^/aho for different values of /3e/ai lo . In 
the small range of energies shown in the figure, which 
corresponds to e ~ Tluj <§; s E , the universal spectrum ap- 
proaches that of the shape-independent results [|J in the 
limit of /3 6 /a ho -40 0. 

Figures |4] and [5] both illustrate the universal p wave 
spectrum at length scale (3®. They also serve to illus- 
trate how different representations of the universal spec- 
trum can serve different purposes in terms of physical 
understanding. In Fig. |4] the universal p wave spectrum 
is plotted as a function of ao//?6 for different values of 
Pe/o-ho- It gives the best illustration that, in the case of 
single channel, the s wave scattering length determines 
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FIG. 3: Universal s wave spectrum at length scale /3e for 
two identical atoms in a symmetric harmonic trap and with 
an asymptotic interaction of the type of — 1/r 6 . The arrow 
points to the s wave scattering length, ao x , beyond which the 
least-bound molecular s state is pushed to positive energies. 



FIG. 4: Universal p wave spectrum at length scale Pe for 
two identical atoms in a symmetric harmonic trap and with 
an asymptotic interaction of the type of —1/r 6 , plotted here 
versus a scaled s wave scattering length, ao/P&. 



not only the s wave spectrum, but also the spectra of 
other partial waves including the p wave [2lj. It also 
determines the p wave scattering length a\ (also called 
scattering volume as it has the dimension of a volume) 
through 



ai/Pl = 



[r(l/4)] 2 a - a 
187T 2do — &o 



(41) 



where a = 2n(3 6 /[T(l/4)} 2 = 0.4779888/3 6 is the mean s 
wave scattering length of Gribakin and Flambaum [33| . 
Figure 0] shows that the p wave spectrum, in the case of 
single channel, is strongly influenced by the atomic inter- 
action only when the s wave scattering length is around 
2do = 0.955978/36, which corresponds to having a p wave 
bound state right at the threshold [13, H{|. For a slightly 
larger than 2a , a\ is large and negative [see Eq. (|4T]) ]. 
the lowest few trap states are strongly affected by a p 
wave shape resonance near the threshold (30| . For ao 
slightly less than 2ao, a± is large and positive, and there 
is a p wave molecular state close to the threshold 3l| . 

In Fig. [51 the universal p wave spectrum is plotted as 
a function of a\/a\ for different values of (3e/ai lo . The 
advantages of this representation are twofold. First, it fa- 
cilitates comparison with the p wave shape-independent 
approximation to be discussed in Sec. HVBl Second, 



while the representation shown in Fig.[4]is applicable only 
in the case of single channel, the representation shown 
in Fig. [5] would apply even in multichannel cases, pro- 
vided the energy dependence of K c due to closed chan- 
nels is negligible in energy range of interest 0, [28l |. 
This is because in representing the p wave spectrum 
as a function of oi, instead of ao, one is only making 
use of the energy-independence of K c , not its angular- 
independence, which takes on different characteristics 
(still related) for multichannel problems [2(1 [28j] . Specif- 
ically, this representation is obtain from solving Eq. (|33|) 
using K c = K c (e — 0,1 — 1), which is related to the ai 
through the following relations [3(| [3l| 



di = —a\ I 1 



1 



*? =1 (e = 0) 



(42) 



where a\ = [r(l/4)] 2 /3|/(367r) is a mean p wave scatter- 
ing length, and 



K?(e = 0) 



Cj ~ K c (e = 0, 1) 
l + ciK c (e = 0,0 



(43) 



where c; = tan(^7r + ^ir). 

Comparing either Fig. 0] or [5] with Fig. pleads to one 
of the more important conclusions of this work. That is, 
unlike the noninteracting particles in a trap for which the 
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FIG. 5: Universal p wave spectrum at length scale f3e for 
two identical atoms in a symmetric harmonic trap and with 
an asymptotic interaction of the type of — 1/r 6 , plotted here 
versus a scaled p wave scattering length (volume) a\/a\ a . The 
arrow points to the p wave scattering length, ai x , beyond 
which the least-bound molecular p state is pushed to positive 
energies. 



lowest p state energy is always greater than that of the 
s states, a strong p wave pairing (a\ ~ a\ Q or greater) 
for interacting particles can lead to a lower energy state 
around the threshold than a s wave pairing. This state is, 
for large and negative a%, ap wave shape resonance stabi- 
lized by the trap, and is a p wave molecular state for large 
and positive a\. Similar statements can also be made for 
other, higher partial waves. Furthermore, we expect the 
same physics to persist in many-atom systems, which will 
be a subject of future investigations. Other characteris- 
tics of the universal spectra are addressed separately in 
subsequent sections. 



B. Limitations of "shape-independent" 
appr oximat ions 



All previous results on two atoms in a symmetric har- 
monic trap can be easily derived as various approxi- 
mations within our theory. In particular, the "shape- 
independent" approximation corresponds to ignoring en- 
ergy dependence of K^ ho \ and taking it to be its value 




FIG. 6: Universal s state spectrum at length scale (3% for 
two identical atoms in a symmetric harmonic trap and with 
an asymptotic interaction of the type of —1/r 6 , plotted to 
illustrate the failure of the shape-independent approximation 
away from the threshold. 



at zero energy, namely, 



3/2)] 2 



Who/2) 



21+1 



lim 



tan Si 
■-p+T 



(44) 



For partial waves with well-defined scattering lengths, 
ai = limfc^o[~ tan<5//fc 2i+1 ], this approximation used in 
Eq. (flT|l leads to the following equation for the energy 
spectrum: 



(-1)' ( I 



1+1/2 



T[l - (e + 1 + 3/2)/2] 



21+1 



2 J r[(-e + Z + 3/2)/2] 

(45) 

for both positive and negative energies. For I = 0, 
Eq. (145)) reproduces the result of Ref. @ , derived using 
a delta- function pseudopotential 

Figure [3] has shown that for the s wave, the shape- 
independent approximation gives a good approxima- 
tion to the universal spectrum under weak confinement 
{(3&/ciho <C 1) and for energies |e| ~ Uui <§; se- The 
shape dependence is more important for strong coupling 
(ao ~ a/i or greater), but especially for energies fur- 
ther away from the threshold. In particular, the shape- 
independent approximation breaks down for all energies 
e ~ se or greater, for which the energy dependence of 
K<- h °) due to the Ion g-range van der Waals interaction 
can no longer be ignored [26|, [3(|. For example, it does 
not give the proper molecular binding energy for small 
positive scattering lengths, and it fails completely to de- 
scribe molecular states of negative scattering length, as 
illustrated in Fig. [S] In reality, these more deeply bound 
molecular states approach those of a free molecule in the 
absence of a trap, with the trapping potential serving 
as a weak perturbation, as to be discussed further in 
Sec. lTVDl 

With the relative success of the "shape-independent" 
approximation for the s wave in the threshold region, it 
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is important to emphasize its severe limitations for any 
partial waves other than the s wave. For atoms with 
—Cq/t & type of van der Waals interaction, the "shape- 
independent" approximation clearly fails for I > 2, for 
which there are no well-defined scattering lengths [30l[34| . 
Even for the p wave, it is applicable only at zero energy, 
or for weak p wave coupling as characterized by a\ ~ 
/3| <C a\ Q or smaller, as illustrated in Fig. |U This failure 
of the "shape-independent" approximation for / > is 
directly related to the failure of the effective range theory 
(ERT) [l5[ in describing the shape resonances close to the 
threshold, and more generally to the failure of ERT in 
describing Feshbach resonances in nonzero partial waves 

Our theory allows for simple generalizations beyond 
the shape-independent approximation. For example, for 
positive energies, Eqs. (fTT|) . (TT2| and (|30|) that determine 
the energy spectrum can be rewritten as 

~^ (/lo) ( e ) = tan^(e). (46) 

Instead of the effective-range expansion for tan 5;, which 
leads to the "shape-independent" approximation, one can 
simply use the corresponding QDT expansion for —Ce/r 6 
type of potential [30]. The results would be applicable 
from zero energy up to e ~ Se for the s wave, and over 
a greater range of energies for higher partial waves. This 
energy range, while much smaller than that described 
by the universal spectrum Eq. (f3"3"|) . already exceeds the 
range of interest in existing experiments [121 ] . 



beyond which the least bound molecular states is pushed 
to positive energy. For a fixed ao (3$, the same equa- 
tion determines the crossing trap frequency 



[r(i/4)r 



(27T) 



■mar. 



(49) 



beyond which the least bound s wave molecular state is 
pushed into positive energy. 
For the p wave, it gives 



aix = [r(1/4)]2 aL = 1.9119552 < o 



(50) 



For a fixed ai 3> it determines the crossing trap 
frequency 



47r 2/3 



[r(i/4)] 4 / 3 ma 2/3 ' 



(51) 



beyond which the least bound p wave molecular state is 
pushed to positive energy. 

All higher branches of states in Figs. [3][5] are trap 
states. All lower branches, which approach those of 
molecular states in the absence of trapping [2l| , as dis- 
cussed in the next section, are molecular states. 



D. Trapping shift of molecular spectrum 



C. Trap states and molecular states 

The states of two atoms in a trap can be classified into 
trap states and molecular states. The former corresponds 
to states that evolve into diatomic continuum as the trap 
is turned off adiabatically (a^ D — > oo). The latter corre- 
sponds to states that evolve into bound molecular states 
in the same limit. Note that the molecular states would 
not have existed in the hard-sphere atomic model Q . 

The molecular state of highest energy correspond to 
the ones in Figs. [3][5] that cross the zero energy. It is 
the least-bound molecular state that gets pushed up in 
energy by the trapping potential. The crossing into pos- 
itive energy can happen either by tuning up the scatter- 
ing length, which has the effect of making the binding 
energy of the molecule sufficiently small in the absence 
of the trap, or by tuning up the trap frequency. 

The crossing points, which correspond to having a 
state right at e = 0, can be easily found through Eq. (|4"5[) . 
which is exact at zero energy for partial waves with well 
defined scattering lengths. For s and p waves, it gives 

2+1/2 



ai x = (-1)' 



r(l/4-i/2) 2l+1 
T(Z/2 + 3/4) ho 



For the s wave with a fixed trap frequency, it gives 

]r(i/4)] 2 



2tt 



-a ho = 2.0920992 a ho , 



(47) 



(48) 



Except for the least-bound state with binding energy 
comparable to or smaller than Tiui, the effects of trap- 
ping on molecular states are generally weak, and can be 
treated perturbatively. The finite range of such states is 
such that the atoms in them would hardly feel the exis- 
tence of a trap in their relative motion. 

The nature of this perturbation is best understood by 
rewriting Eq. ([3"3"|) for e < as 



c(6) 



K C W, 



c(6) 



w, 



:(6) 



(52) 



Comparing this equation to x^ 6 \e s ) = K c , which de- 
termines the molecular spectrum in free space [2lj . the 
effect of trapping is isolated in this formulation to the 
second term in Eq. ([52)) . From Eq. ([37[) and 



Aho) , s -e»i 1 , 1 q - 1/2)0 + 1/2)(Z + 3/2) 

Si ' ' 

it is clear, as expected, that the molecular states with 
binding much greater than fuj are only weakly affected 
by the trapping, and can be treated by solving Eq. (|52[) 
perturbatively. To the lowest order in fte/aho, the energy 
shift due to trapping is given by 



Ae; s = qi(ei s ) [ — 
aho 



(54) 
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where 



m(eis) = (-1)' 



(l-l/2)(J+l/2)(J + 3/2) 



6\e is 



w c(6) w c(fi), 



s=£ls 

(55) 

and ei s — ci/se is the scaled bound state energy of a 
molecule in the absence of the trap. Equation ([54]) means 
that trapping shift is, to the lowest order, proportional 
to (/?6 / 'dho) A , multiplied by a universal function of the 
scaled binding energy that goes to zero in the limit of 
\eu\ -> oo. 



E. Highly excited trap states 



function for two identical atoms in a symmetric harmonic 
trap can indeed be written in the Jastrow form but with 
a correlation function that behaves as 



F(r) 



Cr e 



-3/2 



(59) 



in the limit of large interatomic separations. Since e 
can deviate significantly from 3/2 for strong coupling 
(ao/aho ~ 1 or greater), as illustrated in Fig. [3J this 
result implies that there exists significant long-range cor- 
relation for strongly-interacting particles in a trap. In a 
recent publication [24|, we have shown that such long- 
range correlation exists not only for two particles in a 
trap, but also for N (N > 2) strongly-interacting parti- 
cles in a trap. 



For an highly excited trap state with e ^> 1, the deter- 
mination of energy spectrum can also be further simpli- 
fied. From 

^ o) (e) tan[vr(e - I + l/2)/2)], (56) 
Eq. (pP)]) reduces, for e 3> 1, to 



tan[?r(e - I + l/2)/2] 



tan[<5;(e)] 



(57) 



or (e — / + l/2)7r/2 + 81(e) = jir, where j is an integer. 
This result for I = has also been derived by Bolda et 
al. using a generalized pseudopotential [3] • It could also 
have been derived by using a semiclassical approximation 
for the harmonic part of the potential. For all energies of 
experimental interest, the phase shift Si can be accurately 
described using AQDT parametrization of Ref. [2l[ . 



F. Long-range correlation between atoms in a trap 

In studies of quantum few-body and quantum many- 
body systems, it is often assumed that the wave function 
can be written in a Jastrow form [35j . which is given, 
e.g., for bosons by 



* = 



II 



N 



n 



(58) 



J »<7=1 



Here N is the number of particles, <fi represents an 
independent-particle orbital, and F is the pair correla- 
tion function. It is commonly assumed that F has the 
asymptotic behavior of F(r) — ► 1 at large r, meaning 
that the particles become uncorrelated at large separa- 
tions [351]. 

Our theory here provides a opportunity to check the 
validity of these assumptions, at least for N — 2. To be 
specific, we restrict ourselves here to the lowest s wave 
trap state for which the Jastrow assumptions are usually 
applied. Combining the wave function for the relative 
motion, as given by Eqs. (|5j)- (fTU|) , and that for the center- 
of-mass motion, which is a pure harmonic oscillator in 
its ground state, it is easy to show that the total wave 



V. CONCLUSIONS 

In conclusion, we have presented a QDT for a harmonic 
potential and a two-scale QDT for two identical atoms in 
a symmetric harmonic trap. It is a general theory that 
is applicable to different partial waves, and from deeply 
bound molecular states to highly excited trap states. The 
only approximation in the theory, fte/aho "C 1, can be 
relaxed if needed. The result will still be of the form of 
Eq. (f3"3"]) except that the x^ ho ^ function for 0e/a>ho ^ 1 or 
greater will require a more general solution of the two- 
scale potential described by Eq. (fT9|) . If the scattering 
length parameter used in our formulation is achieved by 
tuning around a Feshbach resonance [3|| [37j , the same 
results would apply, provided it is a broad Feshbach reso- 
nance with a width much greater than the energy scale s e 
associated with the van der Waals potential [H, H^, H(| ■ 

From a more general perspective, the theory we have 
presented here demonstrates the concept of multiscale 
QDT, which has potential implications for a range of 
problems. As examples, we mention two other two-scale 
potentials: 



and 



V(r) = -C 6 /r 6 - C 8 /r 8 



V(r) = -Z/r - d/r 4 



(60) 



(61) 



The former is of interest for a more accurate QDT de- 
scription of long-range molecules over a wider range of 
energies, the latter is of interest for a more systematic 
understanding of the core polarization effect on atomic 
spectra. While such potentials have been well studied 
in numerical calculations, a more systematic approach to 
such two-scale problems may well be worthy of future 
efforts. 
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